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■ Abstract 

In this paper, we use a statistical thermodynamic approach to quantify the em- 
ciency with which clathrates on the surface of Titan trap noble gases. We consider 
different values of the Ar, Kr, Xe, CH4, C2H6 and N2 abundances in the gas phase 
£j I that may be representative of Titan's early atmosphere. We discuss the effect of the 

various parameters that are chosen to represent the interactions between the guest 
species and the ice cage in our calculations. We also discuss the results of varying 
, the size of the clathrate cages. We show that the trapping efficiency of clathrates 

is high enough to significantly decrease the atmospheric concentrations of Xe and, 
to a lesser extent, of Kr, irrespective of the initial gas phase composition, provided 
that these clathrates are abundant enough on the surface of Titan. In contrast, we 
find that Ar is poorly trapped in clathrates and, as a consequence, that the at- 
mospheric abundance of argon should remain almost constant. We conclude that 
the mechanism of trapping noble gases via clathration can explain the deficiency 
in primordial Xe and Kr observed in Titan's atmosphere by Huygens, but that this 
mechanism is not sufficient to explain the deficiency in Ar. 
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1 Introduction 



Saturn's largest satellite, Titan, has a thick atmosphere, primarily consisting of 
nitrogen, with a few percent of methane (Niemann et al., 2005). An unexpected 
feature of this atmosphere is that no primordial noble gases, other than Ar, 
were detected by the Gas Chromatograph Mass Spectrometer (GCMS) aboard 
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the Huygens probe during its descent on January 14, 2005. The detected Ar 
includes primordial 36 Ar (the main isotope) and the radiogenic isotope 40 Ar, 
which is a decay product of 40 K (Niemann et al., 2005). The other primor- 
dial noble gases 38 Ar, Kr and Xe, were not detected by the GCMS instrument, 
yielding upper limits of 10~ 8 for their mole fractions in the gas phase. Further- 
more, the value of the 36 Ar/ 14 N ratio is about six orders of magnitude lower 
than the solar value, indicating that the amount of 36 Ar is surprisingly low 
within Titan's atmosphere (Niemann et al., 2005). These observations seem to 
be at odds with the idea that noble gases are widespread in the bodies of the 
solar system. Indeed, these elements have been measured in situ in the atmo- 
spheres of the Earth, Mars and Venus, as well as in meteorites (Owen et al., 
1992). The abundances of Ar, Kr and Xe were also measured to be oversolar 
by the Galileo probe in the atmosphere of Jupiter (Owen et al., 1999). 

In order to explain the observed deficiency of primordial noble gases in Titan's 
atmosphere, Osegovic & Max (2005) proposed that the noble gases within the 
atmosphere could be trapped in clathrates located on the surface of Titan. 
The authors calculated the composition of clathrates on the surface of Titan 
using the program CSMHYD (developed by Sloan (1998)) and showed that 
such crystalline ice structures may act as a sink for Xe. The facts that the 
code used by Osegovic & Max (2005) is not suitable below about 140 K for 
gas mixtures of interest, and that the authors did not explicitly calculate the 
trapping efficiencies of Ar and Kr in clathrates on the surface of Titan led 
Thomas et al. (2007) to reconsider their results. In particular, Thomas et al. 
(2007) performed more accurate calculations of the trapping of noble gases 
in clathrates using a statistical thermodynamic model based on experimental 
data and on the original work of van der Waals & Platteeuw (1959). On 
this basis, Thomas et al. (2007) showed that Xe and Kr could have been 
progressively absorbed in clathrates located at the surface of Titan during 
its thermal history, in contrast with Ar, which is poorly trapped in these 
crystalline structures. They then concluded that their calculations are only 
partly consistent with the Huygens probe measurements, since the presence of 
clathrates on the surface of Titan cannot explain the primordial Ar deficiency 
in its atmosphere. 

In this paper, we aim to extend the work of Thomas et al. (2007) by con- 
sidering the clathration of a more plausible gas mixture representative of the 
composition of Titan's atmosphere, in which we include Ar, Kr and Xe all 
together. Indeed, Osegovic & Max (2005) and Thomas et al. (2007) both con- 
sidered atmospheric compositions containing only one noble gas at a time. As 
we shall see, the competition between the various species strongly affects the 
efficiency with which they are trapped in clathrates. In addition, because the 
incorporation conditions of guest species in clathrates depend on the structural 
characteristics of the crystalline network and on the intermolecular potentials, 
we examine the influence of the size of cages and of the interaction potential 
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parameters on our calculations. The study of the influence of these parame- 
ters is motivated by two distinct facts. First, it is known that the size of the 
clathrate cages depends on the temperature (Shpakov et al., 1997; Belosludov, 
2002; Takeya et al., 2005), and secondly, several different sets of Kihara pa- 
rameters are available in the literature (Parrish & Prausnitz, 1972; Diaz Peha 
et al., 1982; Jager, 2001). As a result, this study of sensitivity to parameters 
allows us to better quantify the accuracy of the results. 

In Section 2, we describe the statistical model used to calculate the compo- 
sition of the clathrates. This hybrid model is based on the work of van der 
Waals & Platteeuw (1959), and on available experimental data. We also com- 
pare the dissociation pressure of clathrates obtained from our model with that 
obtained from the CSMHYD program (Sloan, 1998). CSMHYD uses a more 
sophisticated and rigorous approach, but is limited to carrying out calculations 
above about 140 K (Sloan, 1998). In Section 3, we investigate the sensitivity 
of our model to various parameters, and also examine the influence of these 
parameters on the predicted clathrate composition. In Section 4, the statistical 
approach developed in Section 2 is used to calculate the relative abundances 
of guests trapped in clathrates that may exist on the surface of Titan. Several 
hypotheses for the abundance of noble gases in the atmosphere of Titan are 
tested. Section 5 is devoted to the summary and discussion of our results. 



2 Theoretical background 

To calculate the relative abundance of guest species incorporated in a clathrate 
from a coexisting gas of specified composition at given temperature and pres- 
sure, we follow the formalism developed by Lunine & Stevenson (1985), which 
is based on the statistical mechanics approach of van der Waals & Platteeuw 
(1959). Such an approach relies on four key assumptions: the host molecules 
contribution to the free energy is independent of the clathrate occupancy 
(this assumption implies in particular that the guest species do not distort 
the cages), the cages are singly occupied, there are no interactions between 
guest species in neighboring cages, and classical statistics is valid, i.e., quan- 
tum effects are negligible (Sloan, 1998). 

In this formalism, the occupancy fraction of a guest species G for a given type 
t (t — small or large) of cage, and for a given type of clathrate structure (I or 
II) can be written as: 



where Cc,t is the Langmuir constant of guest species G in the cage of type 



VG,t = 
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t, and Pg is the partial pressure of guest species G. Note that this assumes 
that the sample behaves as an ideal gas. The partial pressure is given by 
Pq = xq x P, with xq the molar fraction of guest species G in the initial gas 
phase, and P the total pressure. The sum, in the denominator runs over 
all species J which are present in the initial gas phase. 

The Langmuir constants indicate the strength of the interaction between each 
guest species and each type of cage. This interaction can be accurately de- 
scribed, to a first approximation, on the basis of the spherically-averaged Ki- 
hara potential wg{t) between the guest species G and the water molecules 
forming the surrounding cage (McKoy & Sinanoglu, 1963), written as: 




where R c represents the radius of the cavity assumed to be spherical, z is the 
coordination number of the cell and r the distance from the guest molecule 
to the cavity center. The parameters R c and z depend on the structure of the 
clathrate and on the type of the cage (small or large), and are given in Table 
1. The functions 5^(r) are defined as : 

^■5[('-rr-( l+ rD1' (3) 



where ao, 0"g and eg are the Kihara parameters for the interactions between 
guest species and water. The choice of the Kihara parameters for the guest 
species considered in the present study is discussed in the next section. The 
parameters chosen for our calculations are given in Table 2. 

The Langmuir constants are then determined by integrating the Kihara po- 
tential within the cage as 

R 

n 47r f ( w G( r ) \ 2 , ( a\ 

^ = ^r/ exp (-^r) rdr ' (4) 



where T represents the temperature and ks the Boltzmann constant. 

Finally, the relative abundance fa of a guest species G in a clathrate (of struc- 
ture I or II) is defined as the ratio of the average number of guest molecules 
of species G in the clathrate over the average total number of enclathrated 
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molecules : 



fa = 



bsVcs + b e y G ,e 



(5) 



b s E j y.j, s + h E j y.j,e ' 



where the sums in the denominator run over all species present in the system, 
and b s and be are the number of small and large cages per unit cell, respectively. 
Note that the relative abundances of guest species incorporated in a clathrate 
can differ strongly from the composition of the coexisting gas phase because 
each molecular species has a different affinity with the clathrate. 

The calculations are performed at temperature and pressure conditions at 
which the multiple guest clathrates are formed. The corresponding temper- 
ature and pressure values (T = T^jf and P = P£&) can be read from the 
dissociation curve of the multiple guest clathrates. 

In the present study, the dissociation pressure is determined from available 
experimental data and from a combination rule due to Lipenkov & Istomin 
(2001). Thus, the dissociation pressure P^ of a multiple guest clathrate is 
calculated from the dissociation pressure P G 1SS of a pure clathrate of guest 
species G as 



where %g is the molar fraction of species G in the gas phase. 

The dissociation pressure P G 1SS is derived from laboratory measurements and 
follows an Arrhenius law (Miller, 1961): 



where P G 1SS and T are expressed in Pa and K, respectively. The constants A 
and B used in the present study have been fitted to the experimental data 
given by Lunine & Stevenson (1985) and by Sloan (1998) and are listed in 
Thomas et al. (2007). 

The present approach differs from that proposed by Sloan (1998) in the CSMHYD 
program, in which the dissociation pressure of the multiple guest clathrate is 
calculated in an iterative way by requiring that the chemical potential in 
the clathrate phase is equal to that in the gas phase. The determination of 
this equilibrium requires knowledge of the thermodynamics of an empty hy- 
drate, such as the chemical potential, enthalpy and volume difference between 
ice (chosen as a reference state) and the empty hydrate. The experimental 




(6) 



log(J#-)=A+-, 



(7) 
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data available at standard conditions (T = 273.15 K, P = 1 atm) allow the 
CSMHYD program to calculate chemical potentials, and hence dissociation 
curves, as long as the temperature and pressure is not too far from the refer- 
ence point. This method fails to converge at temperatures below 140 K for the 
clathrates considered in this study. Our approach avoids this problem, because 
it uses experimentally determined dissociation curves, which are valid down 
to low temperatures. 

As an illustration, Fig. 1 shows a comparison between the dissociation curve 
obtained from our approach (full line) and that calculated using the CSMHYD 
program (crosses) for a multiple guest clathrate corresponding to an initial gas 
phase composition of 4.9% CH 4 , 0.1% C2H 6 and 95% N 2 . For this compari- 
son, both calculations have been performed with the same set of parameters 
(Kihara parameters and cage geometries given by Sloan, 1998). Both methods 
give very similar results in the range where the CSMHYD program converges. 



3 Sensitivity to parameters 

The present calculations of the relative abundances of guest species incorpo- 
rated in a clathrate depend on the structural characteristics of this clathrate 
(size of the cages for example) and also on the parameters of the Kihara poten- 
tial. It is thus useful to assess the influence of these structural characteristics 
and potential parameters on the calculations, in order to better quantify the 
accuracy of the results. 

3. 1 Discussion of the Kihara parameters 

Papadimitriou et al. (2006) have recently illustrated the sensitivity of clathrate 
equilibrium calculations to Kihara parameters values, in the case of methane 
and propane clathrates. Indeed, by perturbing in the range ± (1%— 10%) the 
a and e Kihara parameters originally given in Sloan (1998), they have demon- 
strated that these parameters have a significant effect on the values of the 
Langmuir constants and on the dissociation curves. It appears thus of funda- 
mental importance to assess the accuracy of the Kihara parameters which are 
used in the studies of clathrates, for example by comparing the theoretical re- 
sults with available experimental data. However, because in the present paper 
we are interested in the calculations of the relative abundances of noble gases 
in clathrates on Titan which are not experimentally available, our choice of 
Kihara parameters has been based on the literature, only. 

For the guest species considered here, i.e., CH 4 , C 2 H 6 , N 2 , Ar, Kr, and Xe, 
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as far as we know, there are only two full sets of Kihara parameters in the 
literature (Parrish & Prausnitz, 1972; Diaz Peha et al., 1982) which have been 
used in the context of clathrate studies. These parameters are unfortunately 
very different (see Table 2). The parameters given by Parrish & Prausnitz 
(1972) have been obtained by comparing calculated chemical potentials based 
on the structural data of the clathrates cages given by von Stackelberg & 
Miiller (1954) with experimental results based on clathrate dissociation pres- 
sure data (Parrish & Prausnitz, 1972). The parameters given by Diaz Peha 
et al. (1982) have been fitted on experimentally measured interaction virial 
coefficients for binary mixtures. They have been recently used by Iro et al. 
(2003) to quantify the trapping by clathrates of gases contained in volatiles 
observed in comets. The corresponding calculations were performed by using 
the clathrate cage parameters given by Sparks et al. (1999). These two sets of 
Kihara parameters (Parrish & Prausnitz, 1972; Diaz Pena et al., 1982) can be 
partly compared to those recently given for CH 4 , C 2 H 6 , N 2 , and Xe, only, in 
the PhD work of M. Jager (2001). 

Table 2 shows that the a and e parameters used by Jager (2001) are quite close 
to those given by Parrish & Prausnitz (1972) and, as a consequence, the rela- 
tive abundances fo (with G = CH 4 , C 2 H 6 , N 2 , or Xe) we have calculated with 
these two sets of parameters are of the same order of magnitude, and behave 
similarly with temperature. By contrast, the relative abundances calculated 
with the parameters given by Diaz Peha et al. (1982) are very different from 
those calculated with the two other sets of Kihara parameters, as expected 
from the conclusions of Papadimitriou et al. (2006). 

Because the potential and structural parameters given by Parrish & Prausnitz 
(1972) (Table 1) have been self-consistently determined on experimentally 
measured clathrates properties, and also because they give results similar to 
those obtained when using Jager's parameters (Jager, 2001) for a reduced set 
of species, we choose here the Parrish & Prausnitz's parameters for the rest 
of our study. 



3.2 Influence of the size of the cages 

In the present paper, we have chosen the potential and structural parameters 
given by Parrish & Prausnitz (1972) (Table 1). However, it has been shown 
that the size of the cages can vary as a function of the temperature (thermal 
expansion or contraction) and also of the size of the guest species (Shpakov 
et al., 1997; Takeya et al., 2005; Belosludov et al., 2002; Hester et al., 2007). 
Indeed, the structural parameter R c increases with temperature and with the 
size of the guest species, whereas it decreases when small guest species are 
encaged. For example, the lattice constant of the methane hydrate is increased 
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by 0.3 % between 83 and 173 K (Takeya et al, 2005). 



In order to quantify the influence of variations of the size of the cages on 
the relative abundances calculated for the multiple guest clathrate considered 
in the present study, we have modified by ±(l%-5%) the R c values given in 
Table 1. These variations are compatible with typical thermal expansion or 
contraction in the temperature range 90-270 K (Shpakov et al., 1997; Takeya 
et al, 2005; Belosludov et al., 2002; Hester et al., 2007). 

Figure 2 shows the evolution with temperature of /at, /kr and /x e (calculated 
from Eq. 5) in a multiple guest clathrate, for both structures I and II, and for 
variations of the size of the cages equal to ± 1% and ± 5%. The calculations 
have been performed for an initial gas phase containing CH 4 , C2H 6 , N 2 , Ar, 
Kr and Xe. The gas phase abundance of CH 4 (4.92 %) has been taken from 
Niemann et al. (2005), whereas the values for C 2 H 6 (0.1 %), Ar (0.1 %), Kr 
(0.1 %) and Xe (0.1 %) are based on our previous study (Thomas et al., 2007). 
The relative abundance of N 2 (94.68 %) has been determined accordingly. 

Our results show that a small variation (± 1%) of the size of the cages has only 
a very small effect on the trapping of noble gases in the corresponding multiple 
guest clathrate, irrespective of the temperature (Fig. 2). In particular, the 
behavior (i.e., increase or decrease) with temperature is not affected by small 
variations of the size of the cages. A similar small effect is also observed for an 
expansion of the cage by 5 %. By contrast, a large contraction of the cages (R c 
decreased by 5 %) leads to strong modifications of the relative abundances fa 
{G =Ar, Kr, Xe) which can vary by several orders of magnitude. 

The evolution of the relative abundances of all guests in the clathrates is 
also given in Fig. 3 as a function of the sizes of the cages, for given pressure 
and temperature. The corresponding calculations have been performed at a 
pressure P = 1.5 bar (i.e., the present atmospheric pressure at the surface of 
Titan), and at the corresponding temperature given by the dissociation curve, 
i.e., T^ s =176 K (see below). Figure 3 clearly shows that, for the given P 
and T conditions, the contractions of the cages have a larger effect than the 
expansions on the relative abundances calculated in clathrates and that these 
effects are also strongly dependent on the interaction parameters between the 
guest species and the cages. As a consequence, the relative abundances in 
clathrates of the noble gases considered in the present paper appear to be 
much more dependent on the size of the cages than those calculated for CH 4 , 
C 2 H 6 and N 2 . This feature can be related to the values of the e parameters 
of the Kihara potential (see Table 2) which are larger for the noble gases 
than for the other species. Moreover, the relative abundances of the smallest 
species (i.e., the noble gases which have the smallest values of the a Kihara 
parameters) are increased when decreasing the sizes of the cages. 
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To summarize, the results given in Figs. 2 and 3 show that thermal variations 
of the cages need to be taken into account if these variations are greater than 
a few percents. Unfortunately, these variations with temperature are often not 
known, except for a small number of specific systems, such as the methane 
clathrate between 83 and 173 K (Takeya et al., 2005) for which the variations 
of the cages have been shown to be very small. In a first approximation, such 
thermal variations will thus be neglected in our calculations of the fa evolution 
on Titan, as a function of the temperature. 



4 Trapping of noble gases by clathrates on Titan 

The statistical approach outlined in Section 2 is used to calculate the relative 
abundances of CH 4 , C2H 6 , N 2 , Ar, Kr and Xe in a multiple guest clathrate 
(structures I and II), as a function of the temperature. As discussed above, 
the interactions between the guests and the surrounding cages have been cal- 
culated by using the Kihara potential with the parameters given by Parrish 
& Prausnitz (1972), and by disregarding the possible influence of the thermal 
variations of the cages. 

The initial gas phase abundance of CH4 (4.92 %) is taken from Niemann et 
al. (2005), whereas three different sets of initial abundances are considered for 
N 2 , C 2 H 6 , Ar, Kr, and Xe in the atmosphere of Titan. The first set of values 
(hereafter case 1) is derived from the atmospheric composition considered by 
Osegovic & Max (2005) and was also used in our previous study (Thomas et 
al., 2007). The second set of values is calculated under the assumption that 
each ratio of noble gas to methane gas in the atmosphere of Titan corresponds 
to the solar abundance (Lodders 2003) with all carbon postulated to be in the 
form of methane (hereafter case 2). The third set of values is calculated under 
the assumption that each noble gas to methane gas phase ratio derives from 
the value calculated by Alibert & Mousis (2007) for planetesimals produced 
in the feeding zone of Saturn and ultimately accreted by the forming Titan 
(hereafter case 3). In each case, the relative abundance of N 2 (approximately 
95 %) has been determined such as £ch 4 + ^Ar + ^Kr + ^Xc + ^n 2 + x c 2 h 6 — 1- 
The initial gas phase abundances for the three cases are given in Table 3. 

Figure 4 shows that the dissociation curves calculated for the multiple guest 
clathrates that form in the three considered atmospheres exhibit a similar 
behavior, although for a given temperature, the dissociation pressure can vary 
by two orders of magnitude from one case to another one (especially at low 
temperatures). However, for a pressure corresponding to the present pressure 
at the surface of Titan (i.e., P — 1.5 bar), the dissociation temperatures given 
by Fig. 4 for the three cases are within a 20 K range, with corresponding 
values T^ s =176, 167, and 185 K, for cases 1, 2 and 3, respectively. These 
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results indicate that the influence of the initial abundances in the gas phase 
is quite weak on the stability of the corresponding multiple guest clathrate. 

Then, we have calculated the variations with temperature of the relative abun- 
dances fa in the multiple guest clathrate considered in the present study, in 
each case. 

Figure 5 shows the corresponding results calculated for a multiple guest clathrate 
of structure I or II. For case 1, this figure shows that the relative abundances of 
Ar, Kr, CH 4 and N 2 decrease when the formation temperature of the clathrate 
decreases, in contrast with the relative abundances of Xe and C2H6 which 
slightly increase when the temperature decreases. This indicates that the effi- 
ciency of the trapping by multiple guest clathrates decreases with temperature 
for Ar and Kr, whereas it slightly increases for Xe. This result differs from that 
obtained in our previous study (Thomas et al., 2007) in which the trapping 
of both Xe and Kr was found to increase when the temperature decreases. 
However, in this previous study, we performed the calculations by assuming 
the presence of only one noble gas in the initial gas phase, the two others being 
excluded. The difference obtained in the present study for a gas phase contain- 
ing the three noble gases, indicates that there is a strong competition between 
the trapping of Xe, Kr, and Ar, when considering that they can be trapped 
all together. Similar conclusions are obtained when considering the two other 
cases, although the trapping of Xe is found to increase much more than in 
case 1 when the temperature decreases. Also, in cases 2 and 3, the trapping 
of Kr appears almost constant in the whole range of temperatures considered 
in the present study. Note that the absolute values of the relative abundances 
are very different for the three cases due to the different compositions of the 
initial gas phase. 

As a consequence, it is much more useful to compare the efficiency of the 
trapping mechanism in each case by calculating the abundance ratios for the 
three considered noble gases Ar, Kr and Xe. Such a ratio is defined as in our 
previous paper (Thomas et al., 2007), i.e., as the ratio between the relative 
abundance f G of a given noble gas in the multiple guest clathrate and its 
initial gas phase abundance x G (see Table 3). The ratios calculated for Xe, 
Kr, and Ar in the three cases considered here are given in Table 4. These ratios 
have been calculated at the particular point on the dissociation curves (Fig. 
4) corresponding to the present atmospheric pressure at the ground level of 
Titan (i.e., P =1.5 bar and T = T^f x s ). 

Table 4 shows that for this particular point of the dissociation curve, the 
relative abundances of Xe and Kr trapped in multiple guest clathrates are 
much higher than in the initial gas phase, irrespective of the initial gas phase 
composition. By contrast, the relative abundance of Ar is similar in gas phase 
and in the multiple guest clathrate. These results indicate that the efficiency 
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of the trapping by clathrate may be high enough to significantly decrease the 
atmospheric concentrations of Xe and, to a lesser extent, of Kr, irrespective 
of the initial gas phase composition, provided that clathrates are abundant 
enough at the surface of Titan. On the contrary, with an abundance ratio close 
to 1, Ar is poorly trapped in clathrates and the Ar atmospheric abundance 
consequently should remain almost constant. 



5 Summary and discussion 

We have extended the work of Thomas et al. (2007) by considering the clathra- 
tion of a gas mixture presumably representing the composition of Titan's at- 
mosphere, where Ar, Kr and Xe are included all together. In this context, 
we have developed a hybrid statistical model derived from the works of van 
der Waals & Platteeuw (1959) and Lipenkov & Istomin (2001), and using 
available experimental data to constrain the clathrates composition. Because 
it has been shown that clathrates equilibrium calculations are very sensitive 
to the guest species - cage interaction potential, we have compared different 
sets of potential parameters existing in the literature. Our calculations were 
performed using the parameters calculated by Parrish & Prausnitz (1972), be- 
cause these parameters form a consistent set for our application to clathrates 
on Titan. We have also assessed the influence of the thermal variations of the 
size of the cages to better quantify the accuracy of the composition predic- 
tion. We show that these variations need to be taken into account if they are 
greater than a few percents. We have then considered several initial gas phase 
compositions, including different sets of noble gases abundances, that may be 
representative of Titan's early atmosphere. We finally show that the trapping 
efficiency of clathrates is high enough to significantly decrease the atmospheric 
concentrations of Xe and, in a lesser extent, of Kr, irrespective of the initial 
gas phase composition, provided that these clathrates are abundant enough 
at the surface of Titan. On the contrary, with an abundance ratio close to 
1, Ar is poorly trapped in clathrates and its atmospheric abundance should 
remain consequently almost constant. Despite the fact that we consider si- 
multaneously three noble gases in the gas phase composition, in contrast with 
Thomas et al. (2007), we obtain the same conclusions: the noble gases trap- 
ping mechanism via clathration can explain the deficiency in primordial Xe 
and Kr observed by Huygens in Titan's atmosphere, but not that in Ar. 

We note that, even if the Visible and Infrared Mapping Spectrometer (VIMS) 
onboard Cassini was able to see the surface unimpeded, the bulk composition 
of Titan's crust is still unknown. Hence, the presence of clathrates on the 
surface of Titan is difficult to quantify. 

Thomas et al. (2007) calculated that the total sink of Xe or Kr in clathrates 
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would represent a layer at the surface of Titan whose equivalent thickness 
would not exceed ~50 cm. The sink of these noble gases in clathrates requires 
the presence of available crystalline water ice on the surface or in the near 
subsurface of Titan. If an open porosity exists within the top few hundreds 
meters in the icy mantle of Titan, by analogy with the terrestrial icy polar 
caps, the amount of available water ice in contact with the atmosphere of Titan 
would thus increase and help the formation of clathrates inside the pores. 
Moreover, in presence of methane clathrate on the surface of Titan, diffusive 
exchange of noble gases with methane might occur in the cavities, thus favoring 
their trapping in clathrates. One must also note that the efficiency of the noble 
gases trapping by clathrates on Titan can be limited by the very slow (and 
poorly known) kinetics at these low temperatures and the availability of water 
ice to clathration. It is then difficult to estimate the timescale needed to remove 
the proposed quantities of noble gases from the atmosphere of Titan. 

To explain the deficiency in Ar in a way consistent with the present results, we 
can invoke the Titan's formation scenario proposed by Alibert & Mousis (2007) 
and Mousis et al. (2007). According to this scenario, the lack of CO in the 
atmosphere of Titan can be explained if Titan was formed from planetesimals 
that have been partially vaporized in the Saturn's subnebula. The vaporization 
temperature in the Saturn's subnebula (~50 K) needed to explain the loss of 
CO in planetesimals ultimately accreted by Titan is also high enough to imply 
the sublimation of Ar and, in a lower extent, that of Kr (see e.g. Fig. 9 of 
Alibert & Mousis, 2007). Indeed, Kr can also be partially trapped in methane 
clathrates formed in the solar nebula (Mousis et al., 2007). On the other hand, 
Xe still remains trapped in planetesimals because its incorporation occurs in 
conditions close to those required for the methane clathration in the nebula 
(Alibert & Mousis, 2007; Mousis et al., 2007). 

It is important to mention that the composition of Titan's today atmosphere 
is almost certainly different from that in past. In particular, the nitrogen 
isotopes in Titan's atmosphere suggest significant mass loss over time. More- 
over, we do not know the outgassing history of methane. In particular, if the 
methane outgassing is recent (Tobie et al., 2006), the lack of infrared opac- 
ity prior to that era must result in freezing out of nitrogen. It is likely that 
the composition, pressure and temperatures in Titan's atmosphere have dif- 
fered significantly in the past, which will affect the composition and formation 
efficiency of clathrates on the surface. 

Finally, we note that there remains the possibility that the noble gas abun- 
dances are telling a story that entirely differs from the scenario we propose, 
in which neither molecular nitrogen nor noble gases were initially accreted in 
clathrates (Atreya et al., 2006; 2007). In this context, Titan would have formed 
from solids produced at such high temperatures that they would have accreted 
nitrogen essentially as ammonia hydrate. Planetesimals formed in such con- 



12 



ditions would be directly impoverished in noble gases since their trapping in 
clathrates require lower temperature and pressure conditions (Atreya et al., 
2006; 2007). 
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Table 1 

Parameters for the cavities. R c is the radius of the cavity (values taken from Parrish 
& Prausnitz, 1972). b represents the number of small (b s ) or large (be) cages per 
unit cell for a given structure of clathrate (I or II) , z is the coordination number in 
a cavity. 



Clathrate structure I II 



Cavity type 


small 


large 


small 


large 


R c (A) 


3.975 


4.300 


3.910 


4.730 


b 


2 


6 


16 


8 


z 


20 


24 


20 


28 
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Table 2 

Three different sets of parameters for the Kihara potential, a is the Lennard-Jones 
diameter, e is the depth of the potential well, and a is the radius of the impenetrable 
core. These parameters come from the papers of (a) Parrish & Prausnitz (1972), 
(b) Diaz Peha et al. (1982) and (c) Jager (2001). 



Ref 


Molecule 


a(A) 


e/MK) 


a(A) 


(a) 


CH 4 


3.2398 


153.17 


0.300 




C2H6 


3.2941 


174.97 


0.400 




N 2 


3.2199 


127.95 


0.350 




Xe 


3.1906 


201.34 


0.280 




Ar 


2.9434 


170.50 


0.184 




Kr 


2.9739 


198.34 


0.230 


(b) 


CH 4 


3.019 


205.66 


0.313 




C2H6 


3.038 


399.07 


0.485 




N 2 


2.728 


145.45 


0.385 




Xe 


3.268 


302.49 


0.307 




Ar 


2.829 


155.30 


0.226 




Kr 


3.094 


212.70 


0.224 


(c) 


CH 4 


3.1514 


154.88 


0.3834 




C2H6 


3.2422 


189.08 


0.5651 




N 2 


3.0224 


127.67 


0.3526 




Xe 


3.3215 


192.95 


0.2357 



Ar 
Kr 
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Table 3 

Three initial gas phase abundances corresponding to the systems considered by Os- 
egovic & Max (2005) (case 1), the solar nebula (case 2), and Saturnian planetesimals 
(case 3). 

Molecule Molar fractions (%) 





case 1 


case 2 


case 3 


Ar 


0.1 


7.1264xl0~ 2 


2.10506 


Kr 


0.1 


3.44xl0~ 5 


1.38xl0" 3 


Xe 


0.1 


3.8xl0" 6 


1.6xl0~ 4 


CH 4 


4.92 


4.92 


4.92 


N 2 


94.68 


95 


92 


C2H6 


0.1 


8.6978 xl0~ 3 


0.9734 
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Table 4 

Abundance ratios of noble gas in clathrates to noble gas in the initial gas phase for 
Ar, Kr and Xe. These ratios are calculated at P = 1.5 bar, and at the corresponding 
temperature on the dissociation curves (see Fig. 4). This temperature is equal to 176 
K for case 1, 167 K for case 2, and 185 K for case 3. 



Case 


Initial molar fraction 


abundance ratio 


abundance ratio 




in gas 


structure I 


structure II 




Ar gas 






1 


O.lxlCT 2 


0.4 


1.5 


2 


7.1264xl(T 4 


0.8 


5 


3 


2.10506xl(T 2 


0.7 


4.6 




Kr gas 






1 


O.lxlO- 2 


7.6 


35.3 


2 


3.44xl0~ 7 


18.3 


143.7 


3 


1.38xl0~ 5 


11.6 


89.5 




Xe gas 






1 


0.1xl0~ 2 


308 


473 


2 


3.8xl0~ 8 


863.4 


2356 


3 


1.6xl0~ 6 


269.2 


948.7 
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1 0" 2 1 1 1 1 1 1 1 1 1 — 

100 120 140 160 180 200 220 240 260 

Temperature (K) 

Fig. 1. Dissociation curves pressure as a function of temperature for a multiple 
guest clathrate corresponding to an initial gas phase composition of 4.9% of CH4, 
0.1% of C2Hg and 95% of N2. The calculations have been performed using either 
the approach discussed in the present paper (full line) or the CSMHYD program 
proposed by Sloan (1998) (crosses). 
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Structure I 



Structure II 




Fig. 2. Relative abundances of Ar, Kr, and Xe in clathrates as a function of tem- 
perature for structures I and II. The solid lines represent the results obtained with 
the parameters of the cages given in Table 1. The dash-dot-dotted and long dashed 
lines correspond to calculations performed with size of the cages increased by re- 
spectively 1 % and 5 %. The medium dashed and dotted lines are results obtained 
with size of the cages decreased respectably by 1 % and 5 % . 



Structure I 



Structure II 




-5 -4 -3 -2 -1 +1 +2 +3 +4 +5 -5 -4 -3 -2 -1 +1 +2 +3 +4 +5 

Variation of cage sizes (in %) Variation of cage sizes (in %) 

Fig. 3. Relative abundances of guests in clathrates as a function of the cage sizes. 
These results have been calculated at P = 1.5 bar, corresponding to a dissociation 
temperature T=176 K. 



22 




100 120 140 160 180 200 

Temperature (K) 

Fig. 4. Dissociation pressures of multiple guest clathrates, as a function of temper- 
ature for the three cases considered in the present study (see text) : case 1 (solid 
line), case 2 (dotted line), and case 3 (dashed line). 



23 



Structure I 



Structure II 




